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We present a fully analytical solution of the dynamics of two strongly-driven atoms resonantly 
coupled to a dissipative cavity field mode. We show that an initial atom-atom entanglement cannot 
be increased. In fact, the atomic Hilbert space divides into two subspaces, one of which is decoher- 
ence free so that the initial atomic entanglement remains available for applications, even in presence 
of a low enough atomic decay rate. In the other subspace a measure of entanglement, decoher- 
ence, and also purity, are described by a similar functional behavior that can be monitored by joint 
atomic measurements. Furthermore, we show the possible generation of Schrodinger-cat-like states 
for the whole system in the transient regime, as well as of entanglement for the cavity field and the 
atom-atom subsystems conditioned by measurements on the complementary subsystem. 

PACS numbers: 42.50.Pq, 03.67.Mn 



I. INTRODUCTION 

The interaction of a two-level system with a quan- 
tized single mode of a harmonic oscillator, named Jaynes- 
Cummings (JC) model []|, is arguably the most fun- 
damental quantum system describing the interaction of 
matter and light. The JC model has found its natural 
playground in the field of cavity quantum electrodynam- 
ics (CQEDVj^ in the microwave [3, Q and in the opti- 
cal regime [J], as well as in other physical systems, like 
trapped ions [5| or Circuit [6| and solid-state [T'l QED. 
Extensions of the JC model to more atoms and more 
modes, externally driven or not, have been developed 
and, presently, we enjoy a vast number of theoretical and 
experimental developments. Unfortunately, a great part 
of them are not easy to handle and most of the inter- 
esting physics has to be extracted from heavy numerical 
solutions and calculations, especially when realistic dis- 
sipative processes are taking into account. 
The advent of quantum information [8|| has been a fresh 
input in the field of CQED Q, reshaping concepts and 
using it for fundamental tests and initial steps in the de- 
manding field of quantum information processing. Here, 
coherence and the generation of entanglement play an 
important role, and in particular the manner in which 
they are affected by the presence of a dissipative envi- 
ronment [lOj. In spite of their relevance, most of the 
attractive quantum models do not enjoy analytical so- 
lutions and the few available ones are not as close to 
realistic conditions as desired. 

In this paper, we consider a system composed by two 
coherently driven two-level atoms trapped inside a cav- 
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ity and coupled to one of its quantized modes. The ex- 
perimental implementation seems to be feasible due to 
the recent advances in deterministic trapping of atoms 
in optical cavities [ll|, [l2|- We show that under full 
resonance conditions and negligible atomic decays the 
system dynamics can be solved analytically also in the 
presence of cavity field dissipation. With the solutions 
at hand we are able to monitor the purity of the cavity 
field and atom-atom subsystems, as well as the entangle- 
ment and decoherence dynamics of the latter. Atomic 
entanglement cannot be generated. However, the quan- 
tum correlations of suitable entangled states can remain 
completely protected in a decoherence-free subspace. In 
case of atomic states outside this subspace, the decays 
of quantum entanglement, coherence and purity are re- 
markably described by the same function. 
In Sec. [IT] we introduce the master equation of the con- 
sidered open system, in Sec. lIIII we present the analytical 
solutions of the system evolution, in Sec. |TV] we study 
the dynamics of entanglement in the atom-atom subsys- 
tem, in Sec. |V]we consider the conditional generation of 
subsystem states, in Scc lVIl we present fully numerical 
results that confirm our analytical developments, and in 
Sec. IVIII we conclude with a summary of our results and 
physical discussions. In the Appendix we discuss the so- 
lution of the master equation presented in Sec. |TT1 



II. OPEN SYSTEM MASTER EQUATION 

We consider a pair of two- level atoms interacting inside 
a cavity with a field mode also coupled to the environ- 
ment (see Fig. [IJ . A coherent external field of frequency 
ud drives both atoms during the interaction with the 
cavity mode of frequency w/ [ij, [iJI ■ The transition fre- 
quency LOa between excited and ground states, \e)j and 




FIG. 1: Two driven two- level atoms interacting with a dissi- 
pative cavity mode. 



\9)j U — li2), respectively, is the same for the two- 
level atoms. This system could be experimentally im- 
plemented with two-level Rydberg atoms in a microwave 
cavity [0] or with three-level atoms reduced effectively to 
two levels interacting with an optical cavity. Relevant ad- 
vances in cooling and trapping atoms in optical cavities 
have been recently achieved [iH,[i2l- Similar dynamics 
could be also implemented in trapped ions interacting 
with a vibrational mode instead of the cavity mode [9| . 

The Hamiltonian which describes the whole system 
unitary dynamics is 
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where il is the Rabi frequency associated with the co- 
herent driving field amplitude, g the atom-cavity mode 
coupling constant (taken equal for both atoms), a (a^) 
the field annihilation (creation) operator, aj = \g)j{e\ 
((tJ = |e)j((7|) the atomic lowering (raising) operator, and 
a^ = \e)j{e\ — \g)j{g\ the inversion operator. 
In the perspective of possible experimental implementa- 
tion of our scheme we must include the effects of cavity 
mode dissipation and the decay of the atomic upper level. 
Therefore, we must solve the following master equation 
(ME) for the statistical density operator p' of the whole 
tripartite system 



/}' = --[n,p'] + Cfp' + CaP (2) 



where 



Cfp' = — [2a/5'a' — a' dp' — p'a'a] 

1 
^aP' = ^E[2'^^'5'a]-a]<7,p'-/)'aja,], (3) 

where k and 7 are the cavity and atomic decay rates, re- 
spectively, and the environment is modeled by a thermal 
bath at zero temperature. 

Changing to the interaction picture the dissipative terms 
remain unchanged and the ME ^ can be rewritten as 



Pi = -t[H-i^Pi] + ^fPi + ^aPi 



(4) 



where Hamiltonian ([T]) has been replaced by the time- 
independent Hamiltonian Tii = Tio + Tii with 






Hi = hg yj ( aja + aja^ 



(5) 



where we introduced the atom-cavity field detuning pa- 
rameter S = uja — ojf, and from now on we consider a 
resonance condition between the atoms and the external 
field {lod = ^a)- We remark that Hamiltonian ([5]) was 
derived by a standard technique, whereas in the case of 
less simple systems, e.g. three-level atoms, more refined 
treatments are necessary to derive a time-independent 
Hamiltonian such as adiabatic elimination or nonlinear 
rotations, as e.g. in ^^. 

The ME ^ can be solved only by numerical techniques, 
as we discuss in Sec lVIl where we present some numeri- 
cal results for the whole system dynamics. On the other 
hand, if we consider negligible atomic decay (7 = 0) it is 
possible to solve Eq. ^ analytically. First of all, we con- 
sider the unitary transformation hi{t) = exp{jr'HQt\ and 
we derive for the density operator p{t) = U{t)pi{t)W {t) 
the following ME: 
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(6) 



where the transformed Hamiltonian can be written as 
2 
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Here, {|+)j, |— )j} [j = 1,2) is a rotated basis connected 
to the standard basis {|e)j, l^)^} via |±)j = MzIELl, jn 
the strong-driving regime for the interaction between the 
atoms and the external coherent field, f2 3> <?, we can use 
the rotating-wave approximation (RWA) obtaining the 
effective Hamiltonian [13l . [l4| 
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Equation ([5]) outlines the presence of Jaynes-Cummings 
((T^d + ajO,') as well as ant i- Jaynes-Cummings (irjd' + 
(Tjd) interaction terms of each coherently driven atom 
with the cavity mode. We remark that in the framework 
of trapped ions a similar Hamiltonian can be found but 
the cavity mode is replaced by the vibrational mode of 
the ions system [Ig] . In the following we will describe the 
solution of the effective ME: 



P = -^[T-(-cs{t),p] 
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III. EXACT SOLUTION AT RESONANCE FOR 
ATOMS PREPARED IN THE GROUND STATE 

In this section we consider the solution of the ME ((9]) in 
the case of exact resonance 6 = 0, cavity field initially in 
the vacuum state, and both atoms prepared in the ground 
state as an example of separable initial state. We leave 
to future treatments the cases of cavity field prepared 
in more general states. The exact solution of the ME 
for any atomic preparation is described in details in the 
Appendix. It is based on the following decomposition for 
the density operator p{t) of the whole system: 

4 4 

m = E (*i^wb')N)oi = E '5«.WN)oi (10) 

where m}U = {I + +), I + ->, I - +), I - -)} is the 
rotated basis of the atomic Hilbert space. Here we report 
only the final expressions for the field operators pij in the 
present case: 
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(11) 



where we introduced the time dependent coherent field 
amplitude 



a{t) =i-^ (l-e-i' 
k \ 



and the function 

/i(i) = exp |-^i + ^ (1 - e-t' 



(12) 



(13) 



We recall that Pji{t) = phit)- We notice the presence 
of single atom-cavity field coherences whose evolution is 
ruled by the function /i (t) J^] , as well as full atom-atom- 
field coherences ruled by /2(i) = ffit), that we shall 
discuss later on. There are also two one-atom coherences, 
and two diagonal terms, which do not evolve in time, 
corresponding to a pure state 

|0)®(| + -) + |-+»^|0)®(|55)-|ee)), (14) 

where we recognize (up to normalization) a Bell atomic 
state |$~). The explanation goes as follows. First of all, 
if we start with atoms prepared in states of the rotated 
basis used in the decomposition of (jlOp . we obtain much 
simpler results due to the structure of Hamiltonian ([5]) 
on resonance and the obvious uninfluence of dissipation 



on the cavity vacuum. Actually, either the field states 
are coherent, |0) (g) | ± ±) i-^ | =F a{t)) (^ \ ± ±), or there 
is no evolution at all for the states |0) (g) | ± =f), showing 
the presence of an invariant subspace for system dynam- 
ics. It is the component of the initial state |0) (X) \gg) 
in that subspace, (fTH) . which does not evolve. We note 
that starting from atoms prepared in the other elements 
of the standard basis we obtain solutions quite analogous 
to (|lip . as can be easily verified from the Appendix. 
In the transient regime {kt ^ 1), the decoherence 

function can be approximated by /i (t) ~ e 2 where 
a{t) = igt, and the whole system is described by a pure 
Schrodinger-cat-like state 
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If we rewrite this result in the standard atomic basis 
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showing the onset of correlations between atomic states 
and cavity field cat-like states. For the terms related to 
the field subsystem we recover expressions analogous to 
those derived in [17[ for a strongly driven micromaser sys- 
tem. Unlike the present system, in that case the atoms 
pump the cavity mode with a Poissonian statistics, inter- 
acting for a very short time such that the cavity dissipa- 
tion is relevant only in the time intervals between atomic 
injections. Furthermore, the cavity field states are con- 
ditioned on atomic measurements. 

At steady state {kt — » cxj) the density operator is mixed 
and given by: 
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with a^^ = 2i^. Interestingly, the steady state has not a 
fully diagonal structure, i.e., it is not completely mixed, 
in agreement with the previous discussion on the time- 
dependent solution. The change from a pure state to a 
mixed one and the degree of mixedness can be we evalu- 
ated by the purity of the whole system /i(i) = Tr[p^{t)]: 
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In Fig. [2] we show ^(i) as a function of the dimensionless 
time kt for different values of the dimensionless coupling 
constant g/k. We see that the purity decays rather fast 
in the strong coupling regime g/k > 1. 




FIG. 2: Purity fi{t) of the whole system density operator 
vs. dimensionless time kt for dimensionless coupling constant 
g/k: (1) 0.05, (2) 0.2, (3) 0.5, (4) 2. 



A. Subsystems dynamics 

Now we consider the time evolution of the atomic and 
cavity field subsystems in the case of both atoms initially 
prepared in the ground state. The cavity field reduced 
density operator Pf{t) — Tra[f){t)] can be derived by trac- 
ing over both atoms: 



Pfit) 



\-a{t)){-a{t)\+m{Q\ + \a{t)){a{t)\ 



(19) 
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Hence the cavity field mean photon number is (N) 
i|a(i)p and its steady state value {N)ss 
results hold for atoms prepared in any state of the stan- 
dard basis. At any time Eq. (fT9)) describes a mixed state, 
whose purity /i/(i) = Trf[pi{t)] is: 
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In Fig. [3] we show the purity /i/ (t) as a function of di- 
mensionless time kt for different values of the ratio g/k, 
showing a better survival of the purity than for the global 
state of Fig. [2] except in the strong coupling regime. 

The reduced atom-atom density operator Pait) can 
be obtained by tracing over the field variables. In the 
rotated basis {| + +)• | H — ), | — h), | )}, we obtain 



Pt{t) = ] 
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(21) 



The presence of six time-independent matrix elements is 
in agreement with the remarks below Ea. (|14p . The purity 
Pait) = Traip^it)] of the bi-atomic subsystem is: 

Pait)^l[3 + Afi{t) + fi{t)]. (22) 
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FIG. 3: Purity iJ.f{t) of the cavity field subsystem vs. dimen- 
sionless time kt for dimensionless coupling constant g/k: (1) 
0.05, (2) 0.2, (3) 0.5, (4) 2. 



Its behavior is quite similar to the one of Fig.[2]for the 
whole system purity. From Eq. (j2ip we can derive the 
single-atom density matrices and evaluate the probability 
to measure one atom in the excited or ground state 



^e,3(i) = i[lT/l(i)] 



(23) 



Quite similar expressions hold for atoms prepared in any 
state of the standard basis. We see that from measure- 
ments of the atomic inversion I{t) — Pg{t) — Peit) we 
can monitor the one-atom decoherence function /i(i) as 
in [l5|. By rewriting the atomic density matrix ([2T|) 
in the standard basis we evaluate the joint probabilities 
Pimit) = {lm\p'a'^\t)\lm) with {l,m}^{e,g}. The cor- 
responding correlation functions at a given time t are: 
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In Figs. [Ifa-c) we show the correlation functions C;,„(i) 
versus dimensionless time kt and coupling constant f ■ 
At steady-state we see that Cggit),Ceeit) i-^ 3/2, that 
is positive atom-atom correlation or bunching, whereas 
Cegit) >— > 1/2, indicating negative correlation or anti- 
bunching. These results generalize those in [l3|, which 
can be derived from Eas. ((24|) in the limit of negligible 
dissipation fct <C 1. We notice that from the joint atomic 
probability Peg (t) = Pgeit) = |[1— /2(i)] one can monitor 
the two-atom decoherence described by /2(t). 



IV. DYNAMICS OF ENTANGLEMENT AND 
DECOHERENCE OF THE ATOMIC SUBSYSTEM 

In Sec. [nl] we discussed the solution of the ME ^ 
for both atoms prepared in the ground state as an ex- 
ample of separable state. In order to describe also ini- 
tially entangled atoms now we consider the general so- 
lution derived in the Appendix for atoms prepared in a 




a) 



b) 



FIG. 4: Correlation functions (a) Cgg, (b) Ccg, (c) Cee vs. 
dimensionless time kt and coupling constant g/k, showing 
atomic bunching and anti-bunching. 



superposition of Bell states |V'a(0)) = X!i 
the coefiicients Ci are normalized as ^^ 

{\Vt)}i=l,..A 

where |<i>*) 



Cilvi), where 
c,P = 1, and 



{|$+), 1$-), |*+), I*-)} is the Bell basis 
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and |*±) 



|e3>±|ge> 



Tracing 



over the atomic variables the solution for the whole den- 



sity operator (Eqs. pH)) and (|A.6|) ) we derive the cavity 
field density operator generalizing Eq. p9|) : 
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For the mean photon number we obtain (N) {t) — (|ci p + 
|c3p)|a;(t)p that is independent of the coefiicients of 
states |$~) and |^~). On the other hand, tracing over 
the field variables we obtain the atomic density matrix. 
We report the solution in the so called magic basis [l8| 
that can be obtained from the Bell basis simply multi- 
plying |j/2) and ji^s) by the imaginary unit: 
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where we omitted the time dependence for brevity. 
To evaluate the entanglement properties of the atomic 
subsystem we consider the entanglement of formation 
epit) [195 defined as 



<^F{t) = log2 



1 + V/1-C2(t), l+v/l-C2(i) 



(27) 



where C{t) is the concurrence that can be evaluated as 
C{t) = max{0, A4(t) - A3(i) - A2(i) - Ai(i)} (A^ are the 
square roots of the eigenvalues of the non hermitian ma- 
trix p^^ {t){p™^ (t))* taken in decreasing order). From 
Eq. (j26p we can derive the probabilities for joint atomic 



measurements in the standard basis 



Pee,ag{t) = - [|ci ^ (1 + /2(t)) + 2|c2 ^ 

+ |c3p(l-/2(i))±4/i(t)i?e(cic;)] 

Peg,aeit) - -^ [|ci ^(1 - ^(i)) + 2|c4p 

+ |c3p(l + /2(i))±4/i(t)i?e(c3c:)]. (28) 

Also, we can derive the density matrix corresponding to 
a single atom and obtain the atomic probabilities gener- 
alizing Eq. (US]): 

PeAt) = ^ [1 ± 2/i(i)i?e(c*C2 + c;ci)] . (29) 

First we consider the case of a superposition of Bell 
states |<i>*) (i.e., Ci = a, C2 = &e*^,C3 = C4 = with a,b 



real numbers). The initial atomic state \gg) (|ee)) can 
be obtained li a — h — ^ and 9 = tt { 9 — 0). We 

find that the concurrence C{t) vanishes for any time and 
every value of g/k, so that it is not possible to entangle 
the atoms. In the case of atoms prepared in a partially 
entangled state we find that the entanglement of forma- 
tion can only decrease during the system evolution as 
shown for example in Fig. [Sfa) in the case a = b = -k^ 
and 9 = it /A. We also see that the progressive loss of 
entanglement is faster for large values of the parameter 
g/k. For atoms prepared in the maximally entangled 
state |$+) (i.e., 9 — Q,a ^ 1,6 = 0) we derive that the 
concurrence simply reduces to ,f2{t) that also describes 
the whole system decoherence. This important point will 
be discussed later. In Fig.ElJb) we show the entanglement 
of formation as a function of dimensionless time kt. For 
atoms prepared in the maximally entangled state I*!"") 
(i.e., 9 — 0,a = 0,6 ~ 1) the concurrence is always max- 
imum {C{t) = 1). In fact, we can see from Eg. ipSj) that 
the atomic density matrix is always the one of the initial 
state, as expected because \^~) is a linear combination 
of the invariant states | H — ) and | — h). 

Starting from a superposition of Bell states |^*) (i.e.. 



kt 10 




kt 10 



b) 



FIG. 5: Entanglement of formation epit) as a function of kt 
and for different values oig/k: 0.1 (1), 0.5 (2), 1 (3), 5 (4). (a) 

atoms prepared in a partially entangled state ' m — — i 

(b) atoms prepared in the Bell states |$^). 

C3 = a, C4 = be^^,Ci = C2 = 0) we obtain analogous re- 
sults. In particular entanglement cannot be generated 
for atoms prepared in states \eg) and l^e), for the state 
|^+) the concurrence is given by f2{t), and the entan- 
glement of state |^~) is preserved during system evolu- 
tion. We find analogous results also for the concurrence 
of atoms prepared in a superposition of |$~) and 1^"*") 
(i.e., C2 — a,C3 = 6e'^, ci = C4 = 0) or in a superposition 
of |$+) and |*~) (i.e., a = a, C4 = 5e'^, C2 = C3 = 0). 
Let us summarize and discuss the main results in 



the case of atoms prepared in entangled states. If we 
consider a superposition of states I*!*") and |^~) (i.e., 
C2 = a, C4 — 6e*^,ci = C3 = 0) we find that the atomic 
density matrix of Eq. (P^ does not evolve. Actually, the 
atomic subspace spanned by I*!*") and |^~) coincides 
with the time-invariant subspace spanned by \ -\ — ) and 
I — h) , so that it remains protected from dissipation dur- 
ing the system evolution. It provides an example of Deco- 
herence Free Subspace (DFS) 20;]. Atomic entanglement 
injected in the system can thus remain available for long 
storage times for applications in quantum information 
processing |21j . 

If we consider a superposition of states 1$+) and |^+) 
(i.e., ci = a, C3 = 6e'®,C2 = C4 = 0) the concurrence is 
given by: 



C{t) = yja'' + 64 _ 2a?-\fl cos(20)/2(t) 



(30) 



and we find the remarkable result that the initial en- 
tanglement is progressively reduced by the decoherence 
bmction h{^)- Note that C[t) = f2{t) for states |$+), 
|^+), and any superposition as al^"*") ± i6|\['+). To un- 
derstand this point let us consider the specific example of 
the initial state |0) (g) 1$^). The evolved density operator 
of the whole system is 

m = l{\- "(*)) (-"(*)i ® I + +)(+ + 1 + 
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In the limit, kt <<^ 1, of short time and/or negligible 



dissipation, where /2(i) — e" 



'2|a(t)|^ 



with a{t) = igt, the 



system evolves into a pure cat-like state where the atoms 
are correlated with coherent states 
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(32) 

For longer times/larger dissipation, the system coherence 
decays as described by the function /2(i). Let us now 
consider the atomic dynamics disregarding the field sub- 
system. The reduced atomic density operator is 



Pa{t) = l[\++){+ + \ + \--){--\ 

+/2(t)(l + +)(--! + ! --)(+ + !)]. 



(33) 



Hence, as the quantum coherence reduces, simultane- 
ously the atoms lose their inseparability, and the state 
becomes maximally mixed (in the relevant subspace) . In 

fact, the atomic purity is given by /ia(i) = — 2^- Hence 
the time evolution of decoherence, concurrence, and pu- 
rity is described by the function /2(i), which can be moni- 
tored via a measurement of joint atomic probabilities (see 



Eq. dMD) 



PeeW = Pgg{t) = -[l + f2{t)] 



P,g{t) = P^^{t)^-[l~f2{t)] 



(34) 



The probability Pee{t) is shown in Fig. [H] as a function 
of kt for different values of the ratio g/k. For kt <^ I, 
when the whole system is in the cat-like state ([5^ . f2{t) 
quadratically decreases as exp{—2g^t^), independent of 
the dissipative rate k. The subsequent behavior is ap- 
proximately an exponential decay whose start and rate 
depend on the atom-cavity field coupling. For g/k > 0.5, 
that is also below the strong coupling regime, we can 
introduce a decoherence and disentanglement rate 



-rD^k\a'"'\ = 2g 



(35) 



that is again independent of (and faster than) k. A 
physical interpretation of this result is that the more 
coupled the two atoms are to the dissipative cavity 
mode, the more effective becomes the decay of both 
the environment-induced decoherence and the initial 
entanglement. For g/k <C 1, that is in a weak coupling 
regime, the exponential decay of coherence and concur- 
rence starts later and its rate, ^d — ^g'^/k, is slower 
than the dissipative rate k. 

We have shown that under full resonance conditions 



separable state can partially entangle by coupling to a 
thermal cavity field |23| . 



CONDITIONAL GENERATION OF STATES 



In this section we seek information about the states of 
one of the two subsystems conditioned by a projective 
measurement on the other one. 

If the system is implemented in the optical domain 
and the cavity field is accessible to measurements, a null 
measurement by a on/off detector implies the genera- 
tion of a maximally entangled atomic Bell state [IJI • In 
the case of atoms prepared in state \gg) (or |ee)), the 
atomic conditioned state will be the Bell state \^~) (see 
(fni) V Analogously, for initial atomic state \eg) or \ge) 
the atomic conditioned state will be \'^~)- 

Now we consider the evolution of the field subsystem 
conditioned by a projective atomic measurement on the 
bare basis {\ee),\eg),\ge),\gg)}. Starting e.g. from the 
initial state |0) \gg), the cavity field will be in the con- 
ditioned states at a given time t (omitted for brevity in 
this section): 




0.375 



FIG. 6: Atoms prepared in the Bell state l^"""): joint atomic 
probability Pee{t) vs. kt for different values of g/k: 0.1 (1), 
0.5 (2), 1 (3), 5 (4). 

it is not possible to generate or increase the initial 
atomic entanglement. This can be explained looking at 
the initial requirements that allow us to write Eq. (|A.l[ ). 
In particular, the strong driving condition il 3> g, the 
resonance condition S — and the choice of negligible 
atomic decays 7 = are necessary to obtain an inde- 
pendent set of equations for the operators pij and to 
exactly solve the system dynamics. In section IVII we 
will show that removing the strong driving condition 
it is possible to slightly entangle the atoms. On the 
other hand, it can be shown [22| that with off-resonant 
atoms-cavity field interaction it is possible to generate 
maximally entangled atomic states also under strong 
driving conditions. Also we recall that, in the resonant 
case and without driving field, two atoms prepared in a 



Pf,{ee,gg) 



2(3 + /2 T 4/i 

4|0>(0| + /2e2|"l^(|a)(-a 



I -a){-a\ + \a){a\ + 
— a) {a 



Pf,ieg,ge) 



T 2/ie— M0)(-aK|0)(a| + 
+ |-a)(0| + |a)(0| 

I — a){—a\ + \a){a\ 



2(1 - /2) 

he'\'^\'(\a){-a\ + \-a){a\\ 



(36) 



Note that in the limit kt -^ 1 the conditioned field state 
is a Schrodinger-cat-like state 



/f,(ee,gg) 



)f,{eg,ge) 



-a)T2\0) + \a) 



^2(e-2|<i|2^4g-|5|2/2 + 3) 
I - (3) - la) 



V2(l-e-2|"l^) 



(37) 



where a(t) — igt. 

The Wigner functions representing the states (|36p in 



phase space at a given time, are: 

2e-2|/3|^ 



in the Lindblad form: 



W 



f,{'^<^,aa) 



iP) = 



7r(3 + /2T4/i) 



2 + 



-2\a\' 



cosh(4|a|/?Ti/?) 



W. 



f,{eg,ge) 



m 



+ /ae^l"!' cos(4|a|i?e/3) 

T 4/iCOsh(2|a|/TO/3)cos(2|a|i?e/?) 

2e~2|/3|2 r 

e"^l"l cosh(4|a|/m/3) 



7r(l - /2) 

/ae^l"!" cos{4\a\ReP) 



(38) 



In Fig. [7] we iUustrate two of the Wigner functions ((38l 




W[pReeJ 



a) 



W[pF,eg] 



b) 

FIG. 7: Wigner function of the cavity field W[/9y_(ee)] (a) and 
W[pf (eg)] (b), for kt = 0.05, and for f = 80 (a)' and f = 40 
(b). 



in tlie transient kt ^ 1 and in the strong coupling regime 
g ^ k. At steady state the Wigner functions are positive, 
corresponding to the states: 



sSS 



1 



PF,{ee,gg) " g 



a^^)(-a^^ 



PF,{eg,ge) 



where a^^ = *^- 



a^^)(a^^|+4|0)(0|] 
^■^)(-a^^'^V|a^^>(a^^l] (39) 



VI. NUMERICAL RESULTS 

To confirm our theoretical analysis as well as to inves- 
tigate system dynamics without the strong driving con- 
dition and including the effect of atomic decay, where 
analytical results are not available, we numerically solve, 
by Monte Carlo Wave Function (MCWF) method ^, 
the ME lU in the resonant case i5 = 0, that we rewrite 



PI = -UrtePl - PlHl) +J2C^PICj (40) 



where the non-Hcrmitian effective Hamiltonian Tie is 



given by 






(41) 



the Hamiltonian Hj is that of Eq.® ioi S ~ , and 

the collapse operators are Ci^2 = V7<5'i,2, C3 = vka. 
We have introduced the scaled time i = gt so that the 
relevant dimensionless system parameters are: 



\l — — k — — 7=—. 
9 9 9 



(42) 



The system dynamics can be simulated by a suitable 
number Nfr of trajectories, i.e. stochastic evolutions of 
the whole system wave function \tpj{i)) [j = 1,2, ...Ntr)- 
Therefore, the statistical operator of the whole system 
can be approximated by averaging over the Ntr trajecto- 
ries,^ i.e., pi{i) ^ ^ E.=; IV',(i))(V',(i)|. 

First we consider negligible atomic decay 7 = to 
confirm the analytical solutions and to evaluate the ef- 
fect of the driving parameter J7. For numerical conve- 
nience we consider the case fc = 1 so that the steady 
state mean photon number assumes small enough val- 
ues. We consider the atoms prepared in the maximally 
entangled states I*!**). We recall that for the state |$+) 
the theoretical mean photon number is {N{t)) = |a(t)p, 
the atomic populations Pe,g{t) — 0.5, the atomic purity 

P-ait) = — Y^' ^^'^ ^^^ entanglement of formation epif) 
is given by (|27|) where the concurrence C{t) coincides 
with /2(i). In Fig. [8] we show e.g. the mean photon 
number and the atomic probability Pg{t). In the strong 
driving limit, fJ = 20, we find an excellent agreement 
with the predicted theoretical behavior. We note that 
we simulated the system dynamics without the RWA ap- 
proximation so that Pg{t) exhibits oscillations due to the 
driving field. We remark that the entanglement of for- 
mation in the case of state |<I>+) evolves almost indepen- 
dently of parameter Q and becomes negligible after times 
gt « 2. In the case of state |$~) epit) decays in a similar 
way for small values of 51, but it remains close to 1 for 
large enough values of the driving parameter as predicted 
in our analysis. 

Finally, we consider the effect of the atomic decay. For 
example, we consider the atoms prepared in the Bell state 
|$~) in the strong driving limit and for the cavity field 
decay rate k — 1. In Fig. [5] we show the mean photon 
number and the entanglement of formation. We see that 
for 7 up to 10~^ the effect of atomic decays is negligible 
and the results of our treatment still apply. For larger 
decay rates the atomic dynamics becomes no more re- 
stricted within the decoherence free subspace. 
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FIG. 8: Effect of driving parameter Q. for negligible atomic 
decay 7 = 0, fc = 1, atoms prepared in the Bell state \^'^) 
((a),(b)) and |$") ((c), (d)), for values of fi: 0.5 (1), 1 (2), 2 
(3), 20 (4). The theoretical functions are the dashed lines. We 
show in (a),(c) the mean photon number {N{t)), and in (b),(d) 
the atomic probability Pg{t) . The number of trajectories is 
Ntr = 500. 
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FIG. 9: Effect of atomic decay 7 in the strong driving con- 
dition Q = 20, atoms prepared in the Bell state |$~) and 
k = 1: 7 = 0.001(1), 0.01 (2), 0.1 (3). (a) Mean photon 
number (iV(t)), (b) Entanglement of formation eF{t). The 
number of trajectories is Ntr ~ 500. 



VII. CONCLUSIONS 

We have described the dynamics of a system where a 
pair of two- level atoms, strongly driven by an external co- 
herent field, are resonantly coupled to a cavity field mode 
that is in contact with an environment. There arc differ- 
ent available or forthcoming routes to the implementation 
of our model. In the microwave regime of cavity QED 
pairs of atoms excited to Rydberg levels cross a high-Q 
superconductive cavity with negligible spontaneous emis- 
sion during the interaction [9] . Difficulties may arise due 
to the (ideal) requirements on atomic simultaneous injec- 
tion, equal velocity, and equal coupling rate to the cavity 
mode. In the optical regime the application of cooling 
and trapping techniques in cavity QED llj allows the 
deterministic loading of single atoms in a high-finesse 
cavity, with accurate position control and trapping times 
of many seconds [12| . In this regime laser- assisted three- 
level atoms can behave as effective two-level atoms with 



negligible spontaneous emission [l5| . On the other hand, 
trapped atomic ions can remain in an optical cavity for 
an indefinite time in a fixed position, where they can cou- 
ple to a single mode without coupling rate fluctuations 
[2a |. These systems are quite promising to our purposes 
and could become almost ideal in case of achievement of 
the strong coupling regime. 

Under full resonance conditions and starting from the 
vacuum state of the cavity field, and for negligible atomic 
decays and thermal fluctuations, we solved exactly the 
system dynamics for any initial preparation of the atom 
pair, thus deriving a number of analytical results on the 
whole system as well as the different subsystems. These 
results are confirmed and extended by numerical simula- 
tions, e.g. including also atomic decays and investigating 
regimes under weaker driving conditions. 

Here we discuss some results mainly concerning the bi- 
atomic subsystem. First of all we find that if atoms are 
prepared in a separable state, no atom-atom entangle- 
ment can be generated, in agreement with previous re- 
sults that showed the classical nature of atom-atom cor- 
relations under purely unitary dynamics [26l |. If atoms 
are initially quantum correlated, their dynamics can be 
quite different as it can be appreciated in the Bell basis. 

If the initial state is any linear combination of the 
states |$~) and |^~), it does not evolve in time. The 
atomic subspace is free from cavity dissipation and de- 
coherence, hence the initial entanglement can be pre- 
served for storage times useful for quantum computa- 
tion/communication purposes [2l| . 

If the atoms are prepared in any superposition of the 
other Bell states !<&"'"), |\1/"'') we find, remarkably, that 
the same function describes the decay of atom-atom con- 
currence, that is an entanglement measure, as well as 
of environment-induced decoherence and purity. Also 
we show that joint atomic measurements allow monitor- 
ing all these fundamental quantities. By generalizing the 
one-atom analysis of Ref. [ISJ we can describe the non- 
trivial decay process. First of all, in a short time/low 
dissipation limit (kt <g; 1) the whole system is in a pure, 
entangled, cat-like state, as e.g. in Eq. ([?!]) starting from 
|0) (g) |<&+). In this initial stage the decay is quadratic 
in time and independent of the cavity dissipative rate k. 
The subsequent behavior can be well approximated by an 
exponential decay. In a weak coupling regime {g/k <C 1) 
the decoherence and disentanglement rate is jd — S>g'^/k. 
In a strong coupling regime {g/k > 1) the rate is jd — 2g, 
i.e., twice the JC coupling frequency. In this regime the 
whole process is independent of the dissipative rate k. It 
is in fact the strong coupling of each strongly driven atom 
to the cavity field, which is in turn entangled with the 
environment, that rules the decay of both entanglement 
and quantum coherence of the bi-atomic subsystem. 
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APPENDIX: SOLUTION OF THE MASTER 
EQUATION 



We illustrate how to solve the master equation ^ in 
the case of exact resonance, ^ = 0. For the density op- 
erator p{t) of the whole system, we introduce the de- 
composition of Eg. fTU)) . Therefore, the ME is equivalent 
to the following set of uncoupled equations for the field 
operators Pij{t) = {i\pi't)\j)- 



Pii = -ig[a^ +a,pii]+ Cfpii 

/5l2,13 = -Jff(a^ + a)pi2,13 + 'C/pi2,l3 

Pi4 = -igia^ + a, P14} + tfpu 

P22,23,33 — •^/P22,23,33 

y024,34 = -«.gP24,34(a^ + o) + CfP2AM 

PA4, =ig[a) + a, P44] H- £/p44 



(A.l) 



Xii = 


igiP 


+r)xii 






- 


l(4-^'w-+|^iO«' 


Xl2,13 = 


-ig 


d d 1 ■ 
d(3 d(3* 2^^ + ^^ 


Xl2,13 


- 


2[^dp /3*5^.+l/3Pjxi2,r3 


Xl4 = 


-2ig 


' d d ' 

8(3 dp* 


X14 


- 


K4-^-|^^|^iO» 


X22,23,33 = 


-K^.V-.^'^'O-- 


X24,34 = 


-^g 


!p-3l^h'"^"'\ 


X24,34 



2 {'^dp ^ d(3* 

X44 = -ig{P + P*)xii 
d 



2 Vdp ^* dp* 



X24,34 



X44. 



(A.2) 



For the initial atomic preparation we consider the general 
pure state |?/'a(0)) — J2i=i^iWi)j where the coefficients 
Ci are normalized as X]i=i l^iP = 1 ^^^d {|^'i)}i=i,..,4 = 
{|$+), |$"), |*+), I*")} is the Bell basis. By this choice 
we can describe atoms in a maximally or partially entan- 
gled state, as well as separable states with both atoms 
in a superposition state |\l/)a(0) = (ai|e)i -|- bi\g)i) (8) 
(a2|e)2 + 621.9)2) where 



aia2 + 6162 aia2 - hh 

ai&2 +bia2 
C3 = ^ C4 



V2 



V2 
ai&2 - bia2 

71 ■ 



(A.3) 



where the brackets [ , ] and braces { , } denote the 
standard commutator and anti-commutator symbols and 
Pj,i{t) — [Piji^)]^- I^ the phase space associated to 
the cavity field we introduce the functions Xij(/3,t) = 
Trf[prj{t)'D{P)]. We note that the functions XiriP^t) 
cannot be interpreted as characteristic functions _27ii for 
the cavity field, because the operators ptj do not exhibit 
all required properties of a density operator. As a conse- 
quence the functions Xy(/3,i) do not fulfill all conditions 
for quantum characteristic functions. Nevertheless, they 
are continuous and square- integrable, which is enough for 



For the cavity field we consider the vacuum state, so that 
for the whole system the initial state is /5(0) = |0)(0| (8) 
|f/'a(0))(f/'a(0)|, corresponding to the functions 



Xy(/3,0) = J2 <^kcn^Wk)H3) exp{|/3|V2}. (A.4) 



Introducing the real and imaginary part of the variable /3 
we can derive from Eq. (|A.2p a system of uncoupled equa- 
tions such that each of them can be solved by applying 
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the method of characteristics [28| . The solutions are: 
Xii,44(/3, t) =l\ci± C3p exp |-^ T a*{t)P ± a{t)P* 

Xl2,13(/3,i) ^ ~-^{ci + C3){C2 T Ci)* fi{t) X 



exp<i-^ + a(t)/3* 



XiiiP,t) =-(ci+C3)(ci-C3)*/2(i)x 

exp l--^ + a{t)*(3 + a{t)(3* 
X22,33(/3,t) = -|c2TC4pexp j ^ 



1 r /jp^ 

X23(/3,t) = -(C2-C4)(c2+C4)*exp j ^ 

X24,34(/3,t) = "2^*^^ TC4)(ci -C3)*/l(i) X 



exp<{-^ + a*(t)/? 



where Q!(t) and /i(i) are defined in Eqs. p2)) . ([13]) and 
/2(0 = /I'lO- From the above expressions we can recog- 
nize that the corresponding field operators Pijit) are: 



/5ii,44(i) = ^Ici ± csl^l T a(i))(Ta(i)| 

/5i2,i3(i) =^(ci+C3)(c2TC4)*-^|-a(t))(0| 
^ e ^ 

/5i4(t) =^(ci+C3)(ci-C3)*^3§^|-«W)(«WI 

P22,33(i) =i|c2TC4p|0)(0| 

P23W =^(C2-C4)(C2+C4)*|0)(0| 
P24,34(i) = -7^{C2 T C4)(C1 - C3)* }^^^ma{t)\. 

^ e 2 

(A.6) 

Hence, by remembering that pji — p^-, we can recon- 
(A.5) struct the whole system density operator p{t) of Eq. ^ 
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